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ABSTRACT 


This report summarizes the work performed during the first six 
months of a NASA-Lewis sponsored research program (Contract Number 
nag 3-67). This research program is concerned with the calculation of the 
far field radiation patterns surrounding various practical jet engine inlet 
configurations under different excitation conditions. This initial phase of the 
research program is primarily concerned with the upgrading of existing 
computer codes. These computer codes are based on a special integral 
representation of the external solutions of the Helmholtz equation and they 
are capable of calculating the sound field radiated from general 
axisymmetric bodies with complex boundary conditions up to non- 
dimensional wave numbers of ten (based on radius) with less than ten percent 
error in the acoustic quantities of interest. 

During the past six months the computer codes have been refined and 
expanded so that they are now more efficient computationally by a factor of 
about three and they are now capable of producing accurate results up to 
non-dimensional wave numbers of twenty. Computer programs have also 
been developed to help generate accurate geometrical representations of the 
inlets to be investigated. This data is required as input for the computer 
programs which calculate the sound fields. This new ’’geometry genera ting"i 
computer program considerably reduces the time required to generate the 
input data which was one of the most time consuming steps in the process. 

The results of sample runs using the NASA-Lewis QCSEE inlet are 
presented and comparison of run times and accuracy are made between the 
old and upgraded computer codes. The overall accuracy of the computations 
is determined by comparison of the °f the computations with simple 

source solutions. 
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INTRODUCTION 


This research program consists of two major tasks. The first task is 
concerned with the upgrading and testing of the computer codes and the 
second is concerned with the calculation of the radiation patterns for 
specific inlet configurations. This report describes the progress made to date 
on the completion of the first task. 

The computer programs previously developed are capable of 
predicting the sound field radiated by various geometrical configurations. 
The theoretical basis for these programs is presented in Refs. 1>3 (See 
Appendix A for Ref. 3). These publications describe the development of a 
unique integral formulation of the acoustic radiation problem and its 
application in the calculation of 2-D, 3-D and cyiindricaily symmetric 
radiation patterns. In the present investigation the cyiindricaily symmetric 

3 

form of the equations is employed as the inlet geometries we are concerned 
with in this study are all axisymmetric. 

The computer codes developed for the calculation of radiation 
patterns about an axisymmetric body were initially developed for 
calculations at relatively low non-dimensional wave numbers (i.e., less than 
10) and thus were not as efficient as they could be. To calculate results for 
higher wave numbers (i.e., Up to 20) the computer codes had to be upgraded 
so that more calculation points could be handled on the bodiesj and,, to keep 
the calculation times reasonable the computer codes had to be made more 
efficient. Also, in using these "olcf* corhputer codes considerable time was 
required lor the specification of the body geometry and boundary conditions 
(i.e., the input data). To expedite the input process computer codes had to be 
developed to generate the input data in a form which could be directly input 
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into the *'new" computer codes. | 

■ j 

Once the computer codes were upgraded they were tested for speed ! 

and accuracy. To test for speed the upgraded programs were simply | 

compared to a bench mark set of cbmputer runs done with thd old computer 
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codes for the ”QCSEE” Inlet previously studied for the Air Force . The 
accuracy of the computer codes was checked by comparison of the 
calculated solutions to those generated by simple sources. The method of 
using a simple source to generate and check solutions for complex 
geometries is detailed in Ref. 3(See Appendix A.). 

The particular form of the equations used in this study is termed 
cylindrically symmetric as the solutions describing the radiation patterns 
from the axisymmetric inlets may have tangential dependence. In this case, 
the radiation patterns may be calculated separately for any tangential mode; 
and, since each tangential mode must be solved for separately there is no 
coupling between them. Also, since the equations are linear, as many 
tangential mode, solutions as required may be calculated and these modal 
solutions may be summed. For each tangential mode all radial modes (both 
cut-on and cut-off) are automatically considered and thus the theory does 
allow coupling between the various radial modes present for each tangential 
mode. That is, if a particular radial mode is present at the driver (i.e.,the 
fan) plane but it i$ above cut-off for the duct the theory shows that it dies 
out exponentially as it should. 

In this study tangential mode numbers, M, of up to 13 are considered 
with corresponding non-dimensional wave numbers of up to 20. Since the 
cut-off frequency for the first radial mode of the 13th tangential mode is 
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approximately ka = 14.92837 (where k is the wave number and a is the 
characteristic body dimension — typically the duct radius at the driver 
plane) the cut-off ratio is of the order of 1.3. The limiting factor is the non- 
dimensionall wave number so that higher cut-off ratios may be achieved for 
lower tangential mode numbers. In this connection it should be pointed out 
that during the first half of this contract period computer programs have 
also been developed to calculate the non-dimensional cut-off wave numbers 
for any tangential and radial mode and to calculate the mode shapes for 
these conditions. 


EFFICIENCY AND ACCURACY OF THE COMPUTER CODES 


Tile original computer codes were used to produce bench mark runs 
(or comparison purposes. One of these runs (which is representative) was 
made for the case where there were 105 points on the surface of the QCSEE 
inlet (See Appendix A for geometry.) at a non-dimensional wave number of 
ka = 10.0 and a tangential mode number of M s 10. The acoustic velocity for 
a simple source was specified on the surface of the body and the acoustic 
potential was calculated. The calculation of the distribution of the acoustic 
potential on the surface of the body required 867 sec. and the average 
absolute errors in its modulus and phase were found to be and 33.0% 
respectively. With this data available a separate program is then used to 
calculate the acoustic potential at points in the field surrounding the body. 
Since this program involves a straight forward summation; not much 
optimization could be done on it. Thus, most of the optimization work was 
directed towards the program which calculates the acoustic quantities on 
the surface of the body. This surface program contains an NxN matrix 
(where N is the number of points on the body) that must be filled and solved 
to generate the surface quantities. Therefore, the computing time for the 
surface quantities go up roughly as the square of the number of points on the 
body which can become very significant for large N . 

Much of the work in this area was directed simply towards generating 
a more efficient computer code, and significant gains were made. Also, 
different numerical integration methods were tested along with different 
matrix solving schemes; however, these were not as fruitful. For example. 


different integration schemes such as Filon integration^ were tried in the 
circumfrential, 6 , direction* Compared to the presently used Gauss- 

Legendre integration scheme they were found to be no better and in most 
cases worse. Also, different matrix solvers such as Gauss>Seidel iteration 
with both over and under relaxation were tested and it was found that in 
order to maintain the accuracy obtained with the presently used Gauss- 
Jordan matrix reduction scheme too may iteration steps were required. 

Some very minor changes to the computer codes such as adding the 
normals to the body at all the integration points within an integration sub- 
interval, yielded dramatically better results. It was also found that 
specifying the distance in the Gauss-Legendre integration formula in the 
6 direction in reverse order (i.e., starting at 180*^ Instead of 0^) yielded 
better results. Neither of these charges increased computing time 
significantly but they were found to be responsible for much better 
accuracy. 

New computer runs were then made for the same cases as the original 
bench mark set. For the represen titive runs mentioned at the begining of 
this section (i.e. ka s 10,0, M = 10, 105 points on the surface of the QCSEE 
inlet with the normal acoustic velocity specified everywhere) the run time 
required for the surface solution of the acoustic potential decreased by 
almost a factor of 3 from 867 to 309 sec. Also, the average absolute percent 
errors in the modulus and phase of the computed acoustic potential on the 
surface of the inlet decreased from 14.996 and 33.0% to 9.<f% and 10.3% 
respectively. ^ ^ ^ ^ ^ ^ ^ ^ ^ 
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UPGRADING OF THE COMPUTER CODES 

The second effort during the reporting period was concerned with the 
upgrading of the computer codes so that better accuracy could be obtained 
at higher non-dimensional wave numbers. This was done by the addition of 
more points on the surface of the body. It was expected that this effect 
would be linear in that twice the number points on the body should halve the 
error in the potential calculated on the body. 

Another set of computer runs were then performed with 156 points on 
the surface of the QCSEE inlet. The increase in the required computer time 
should fall between 1/3 N^, the operation count for the matrix solution via 
Gauss-3ordan reduction, and N , the number of terms that must be 
calculated to fill the matrix. This being the case the run time should fall 
between 681 sec. and 1013 sec.. For the representative case of M = 10 and 
ka = 10.0 the run time was in fact found to be 791 sec. which is less than the 
original bench mark runs of 867 sec^ for the same case with only 105 points. 

The average absolute percent errors in the modulus and phase of the 
potential on the surface of the body were found to be 8.7% and 93.5% 
respectively, In this connection one should note that although the average 
absolute percent error in the modulus of the acoustic potential decreased it 
did not decrease as much as expected and the error in the phase actually 
went up 9 times. 

This anomoly was found to be due to the way in which the error was 
calculated; that is, the percent error was very small where the values of the 
potential were high but with the same magnitude error the percent errors 
were very high when the potential was small. Upon studying the results of 


various runs it was found that the magnitude of the error remained almost 
constant over the entire body and thus a better measure of the error would 
be the normalized average absolute error in percent. It should be pointed out 
that this will be a departure from the way in which the errors were 
calculated in the past (i.e.,in the proposal for this research) and therefore 
previous error estimates can not be compared directly with these. It is also 
felt by the authors that this is a more reasonable way to calculate the error 
as the high tangential mode numbers whose characteristics are of interest in 
this study have solutions that tend to vary over many orders of magnitude. 
As an example of this, the solution for the acoustic potential on the surface 
of the inlet for the simple source case where M = 10 and ka = 10,0 varies 
from zero at the cent?}| iine of tne inlet to (10^^) near the fan plane on 
the centerbody. The formula for the normalized average absolute error in 
percent is given by: 

E = 1 (tPEXaCT “ 'PcALcM '^n) 

'i“l 

( 1 ) 

These errors, for the modulus and phase of the acoustic potential 
calculated on the surface of the inlet, were then recalculated for the 
representative bench mark run (i.e., M = 10, ka = 10.0, 105 points on the 
surface of the body) and were found to be 0.139% and 8.75% respectively. 
These errors were also calculated for the results generated by the improved 
computer code again using 105 points on the surface of the body and they 
were found to be 0.0999% and 3.10%, respectively. Finally, these same error 
estimates were made for the case where 156 points were used on the 
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surface of the body and they were found to be 0.0^79% and 2.99%. As can be 
seen these error estimates show that the error is decreased by about a third 
as is expected when the number of points is increased by about a half. These 
results are presented in compact form in Table I. 

A run was then made for a case where the cut-off ratio was greater 
than one. The case chosen was for a tangential mode number of 5 and a cut- 
off ratio of 2. This yields a non-dimensional wave number of ka s 12.83. To 
obtain an exact solution for comparison purposes, both on the surface and in 
the field, the artifice of a simple source within the body was again 
employed. Having Calculated the distribution of the acoustic potential on the 
surface of the QCSEE inlet using the normal acoustic velocity distribution is 
input the normalized average absolute errors of the modulus and phase of 
the potential were found to be 0.^f92% and ^.14% respectively. 

The computer program that calculates the potential and the acoustic 
velocity in the field was then used to calculate these acoustic quantities at 
38 points in the field surrounding the inlet. These points were placed every 9 
degrees from the centerline to 90° from the centerline on two quarter 
circles centered at the entrance plane of the inlet at distances (radii) of 25a 
and 100a, where a is the non-dimensionalizing body length (i.e., the inside 
radius of the inlet at the fan plane). Results were then calculated lor two 
different distributions of the acoustic potential on the surface of the inlet. 
The first distribution used was the exact one that had been calculated for 
Comparison purposes and the second one was the distribution calculated 
using the surface code. The exact normal acoustic velocity distribution Was 
used in both cases as it was used as the boundary condition for the surface 
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3 olution. The respective computing times lor 38 points in the field with 138 
points on the suflace of the body are 138 seCt and 148 sec.. The difference in 
computing times is due to the fact that the COC 70/74 computer that was 
used has a time sharing system and thus run times vary slightly as the load 
(number of jobs) on the system changes. 

Using the exact distribution of the acoustic potentiai on the surface 
of the inlet yielded results with normalized average absiute errors in the 
modulus and phase of 0.229% and 30.1%. These results are presented in Figs. 
1 and 2 in terms of both SPL and PWL in dB along with the exact solution in 
the field. The equations used for SPL and PWL in dB are presented below in 
terms of the acoustic potentiai 9 and the outward normal acoustic velocity 

an • 

SPL(dB) « 20 logjjj 

( 2 ) 


PWL(dB) - 10 logj^j k<lcp U > + 143.3 

(3) 

whereo means time average. 

Using the calculated surface solution, the errors were 0.385% and 
36.3% respectively and these results are presented In Figs. 3 and 4. Since the 
error in the phase is so high (compared to the error in the phase in the 
surface solution) a run was made with more points on the surface of the 
body using the exact solution on the surface to see if the results would get 
better as they should. For the case where 230 points were used on the 
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INLET 



Figure 1. Exact Solution Employed on the Surface at 156 
Points, Cut-Off Ratio = 2. 
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surface of the inlet with the exact solution specified for both the acoustic 
potential and normal acoustic velocity the errors in the acoustic potential 
at 38 points in the field were found to be 0,0127% and 1.18% for the modulus 
and phase, respectively^ These results are presented in Figs, 3 and 6 . This 
shows that the calculated field solutions do converge to the proper values as 
more points are specified on the surface of the body. 

Another set of runs, surface soiution and field solution, were then 
done for a tangential mode number of 2 and a cut-off ratio of 3 so that the 
non-dimensional wave number is ka = 15.27. The same number of points were 
used in this case as in the previous one (i.e,, 136 points on the body and 38 
points in the field) and again the artifice of a simple source within the inlet 
was employed. 

For the surface solution the normaiized average absolute errors of 
the modulus and phase of the acoustic potential was found to be 1.112% and 
5.33%, respectively. Since this run is for a relatively low tangential mode 
number M = 2 the modulus of the acoustic potential on the surface of the 
body doesn't change much (i.e., it only changes ^ orders Of magnitude) so 
that the average absolute percent errors are a relatively good measure of 
the error also. For comparison purposes they were computed and found to be 
9.87% and 13.7% for the modulus and phase, respectively. 

The acoustic potential and the outward normal acoustic velocity were 
then calculated in the field on the same two quarter circles as before. The 
normalized average absolute percent errors are found to be i.27% and 8.59% 
respectively for the modulus and phase of the acoustic potential. These 
results are plotted in Figs. 7 and 8 in terms of SPL and PWL in dB as before. 
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Figures. Exact Solution Employed on the Surface at 230 
Points, Cut-Off Ratio = 2. 
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Figure 6. Exact Solution Employed on the Surface at 230 
Points, Cut-Off Ratio = 2. 








QCSEE INLE 


(tS 

(0 o 

ID O 
CV2 


^ O 
c (d 
D X 
X (y. 

Id 

D O 
ID O 




< \- 


o ocoo 
o 


CO O m 0) O !H -St* CO N CO O O -h cm C\3 c\j c\j 

05 O O O O ri ^ ri T-( ri CM CM CM CM C\i CM CM 


>-J 5 OrfCM CO CDO-* CM lO t'- CO 05^ CM CM CO ^ ^ rj* 

0. ID o»-'*-'-iCMCMCMCMCMCMCMC0C0C0C0C0C0C0 

W r-< r-t r-( 1 -x rH rH r-4 rH rl rH 


O ID O D O D O D O D O D OlD O DO to O 
ri CM CM CO CO 'i' ^ D D 50 CD r- N CD CO 05 


Figure 7. 156 Pointson the Surface, Cut-Off Ratio 



QCSEE IN LEI 




Al$o, a summary of the results for the two cases where 1 % points were used 
on the surface of the body are presented in Table II, 

The final two runs were made for a case that will be run for the 
NASA Lewis 3T1SD bellmouth inlet. This particular case was chosen as it 
requires the highest non-dimensional wave number and therefore should 
represent the worst case (i.e.f the case with the most error). For this case a 
tangential mode number of M s 13 and a cut-off ratio of 1.30 were used 
which translates into a non-dimensional wave number of ka = 19.<H. On both 
of the following runs 181 points were specified on the surface of the QCSEE 
inlet. 


For the first run of the set the artifice of a simple source was again 
used so that an exact solution Would be known and the error could be 
calculated. The normalized average absolute error in the modulus and phase 
of the acoustic potential were found to be 0.105% and 24,0% on the surface 
and 0. 1 1 1% and 38.2% in the field. The results in SPL and PWL in dB are 
presented in Figs. 9 and 10. In evaluating these results, it should be pointed 
out that considering the range of dB values covered by the solution, the 
errors in dB are relatively small. 


For the final run, the simple source was not used; instead, the 

acoustic velocity was specified for the 13th tangential mode at the fan plane 

and the rest of the body was specified as having a hard wall (i.e., - o ), 

on 

First, the surface distrubution of the acoustic potential was calculated (1131 
sec.) and then the values of the acoustic potential and outward normal 
acoustic velocity were calculated at 38 points in the field as before (142 
sec.). Finally the results were plotted as before (See Figs. 1 1 and 12.) except 
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Source Distribution at the Fan Planej 
tints on the Syrjface, Cut-Off Ratio = 













that for this case no exact solution was available. It should be noted here 
that a computer program has been developed to calculate the normal 
acoustic velocity distribution at the fan plane for any mode (tangential and 


SUMMARY 


During the first half of this contract period the basic computer 
programs have been made much more efficient and accurate. Furthermore, 
they have been upgraded so that they can handle many more points on the 
surface of the body thus significantly extending their range of applicability. 
Also, new computer programs have been developed to help generate the 
required input data for these programs. This data consists of the body 
geometry and the input modal distribution at the fan plane. Another program 
has been developed which can accurately calculate the modal cut-off wave 
numbers for any tangential and radial mode. Finally, a computer program 
was developed that is capable of plotting the calculated data in the field in 
terms of SPL and PWL in dB (See Figs. 1-12.). 

Various results are calculated for the NASA Lewis QCSEE inlet as 
part of the test program for the ’’new", upgraded computer codes. Tangential 
mode numbers of up to 13 were used with corresponding non-dimensional 
wave numbers over 19. The artifice of a simple source solution was 
employed for most of the computer test cases so that error estimates could 
be calculated. Large improvements in both accuracy and computing time 
required for the "new" computer codes when compared to the "old" 
computer codes were found. 
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Prediction of the sound field radiated from axisymmetric 
surfaces 
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A general analytical melluHl for determining the radiated sound fields from axisymmetric surfaces of 
arbitrary cross tSectioiV with general boundary conditions Is developed* The method Is based on an integral 
representation for the external solutions of the Helmholtz equaliom An integral equation is developed 
governing the surface potential distribution which gives unique solutions at all wavenumbers^ The 
axisymmetric formulation of the problem reduces its solution to the numerical evaluation of line integrals 
by Gaussian quadrature. 1'he applicability of the solution approach for both a sphere and finite cylinder is 
demontrated by comparing the numerical results with exact analytical solutions for both discontinuous and 
continuous boundary conditions. The method Is then applied to a jet>engine>inlet configuration and the 
computed results are in good agreement with exact values, 

VACS numberst 43.20.Rz, 43.20,Tb 


equation in the surrounding infinite three-dimensional 
domain* 


INTRODUCTION 

To reduce the noise radiated to the community from 
turbofan inlets^ the effects of sound suppression mate- 
rial in the Inlet and the spatial distribution of the sound 
source on the radiated sound levels and patterns must be 
determined, Arialytical techniques for predicting these 
effects must be capable of dealing with general axisym- 
metric geometries and complicated boundary conditions 
which are encountered in multiply lined inlets. For in- 
stance^ in a typical inlet the compressor-fan combina- 
tion represents a noise source with a nonuniform spatial 
excitation p«ittern. Thus^ the analytical method should 
be capable pf taking into account sound sources of gen- 
eral spatial dislrlbution , AlS 0 | inlets may contain mul- 
tiple acoustic liners to reduce the radiated sound power/ 
and admittance boundary conditions are commonly used 
to account for the absorption characteristics of the liner. 
Therefore, the analytical method must be capable of 
dealing with spatially varying surface admittances^ 
Finally, the method should be capable of predicting the 
characteri$tics of the radiated sound field in an infinite 
domain. Keeping these requirements in mind, the work 
presented in this paper describes the results of an in- 
vestigation which has been concerned with the analytical 
determination of radiated sound f^Mds from axisymmet-^ 
ric surfaces of arbitrary cross section and with general 
boundary conditions . 

The method used in this Investigation is based on an 
integral form of the solutions of the Helmholtz equa- 
tion.**^ With this formulation the acoustic potential 
anywhere external to the surface can be found once the 
distribution on the surface is known. Thus, to deter- 
mine the radiated sound field the problem reduces to 
the determination of the distribution of the acoustic po- 
tential on the two-dimensional surface of the geometry 
under consideration instead of solving the Helmholtz 


It has been previously shown*"* that when applied to 
exterior sound radiation problems the classical tech- 
niques fail to produce unique solutions at frequencies 
corresponding to certain interior eigenvalues of the ge- 
ometries under consideration. Unless special precau- 
tions are taken, straightforward numerical solution of 
the integral equation produces large errors at frequen- 
cies close to these eigenvalues. For the general geom- 
etries of interest in this study, these eigenfrequencies 
are not known a priori. Therefore, the frequencies 
about which large numerical errors can occur cannot be 
easily avoided . A critical review of available analytical 
techniques for avoiding these errors is provided by Bur- 
ton in Ref. 1. In a search for an appropriate technique 
for use in the present study of inlets, the authors pro- 
grammed each of these methods for a sphere and ob- 
tained numerical results for the surface and radiated 
sound fields. This study showed that the metliod of Bur- 
ton and Miller* was the most straightforward to imple- 
ment. However, an interpretation of a strongly singular 
integral, given in the analysis in Ref. 5 by Meyer ei ql, , 
was necessary for the equations to be amenable to nu- 
merical solution. Basically the method proposed by 
Burton and Miller involves a reformulation of the inte- 
gral equation for the acoustic potential, and the solu- 
tions obtained are valid at all frequencies. It also yields 
the most consistently accurate results for a given number 
of points at which the acoustic potential is numerically 
evaluated on the surface . Tlierefore, the method based 
on the analysis in Ref. 5 has been chosen for this inves- 
tigation. 

The resulting integral equation for the surface acous- 
tic potential is solved numerically and, for axisymmet- 
ric geometries, the equation reduces to the evaliiation 


63 T Acoust. Soc, Am. 65{3), Mar. 1979 0001 *4966/79/030631 *08$00.80 © 1979 Acoustical Society of America 631 


of lino lntcgt*tlj| . Thus > Iho axlsymmctrlc case can be 
reduced id an ci:iulvalcnt one-dlmenslonal problem* Hav- 
ing dlscreUaed the Integral equation, the resulUng sys-^ 
tern of algebraic equations Is solved using complex 
Gauss-Jorcian oltmlnatton. Since the coefficient matrix 
involyes the free^space Green's function, which becomes 
singular as two iKjlnts on the surface approach one an- 
other. luunerical teehmques are proseivted which can 
deal with these singularities and yield accurate results . 
Gaussian integration is used to increase the accuracy of 
the solution without significant penalties in computer 
storage and lime requirements, T(ie applicability of 
the integral formulation and the accuracy of the numeri- 
cal techniques are demonstrated by computing the sur- 
face and farfleld distribution of the acoustic potential 
on both a sphere and a finite cylinder* The numerical 
results are compared with known exact solutions gen- 
erated by the separation of variables technique. Sur- 
faces with spatially varying forcing functions and ad- 
mittances are considered, for different tangential modes, 
to evaluate the capability of the integral approach to 
handle boundary conditions of a general nature. With 
the sphere* agreement between computed and exact re- 
sults is to three signlficanl flpreSt For the cylinder 
agreement 1$ to two significant figures* The effect on 
the accuracy of discontinuous boundary conditions In- 
volving nonzero admittances over the surface and of the 
corners eiveountered In the cylindrical configuration are 
also prosentedi Finally, the numerical results Tor an 
Inlel configuration are compared with exact solutions and 
agreement is to within 10%. 



yiG^ 1, Geomotricul protxiiiios of the general ncousUo rndia- 
tton problem. 


lion of the Helmholtz equation and Is taken to be the Tree- 
space Green's function for a point source® defined as 

C{P,Q)^exp[tkr(P,Q)]/r{P,Q) . (4) 


I. THEORETICAL CONSIDERATIONS 

In this section the general three-dimensional integral 
representation of the solutions of the Helmholtz equa- 
tion is developed for application to radiation problems . 
This particular formulation yields unique solutions at 
all frequencies and does not have strong singularities 
which are difficult to handle numerically. The general 
Integral equation is then specialized for axisymmetrlc 
geometries . A more detailed development is presented 
In Ref. 5* 


A. General formulation 

Beginning with the three-dimensional Helmholtz equa- 
tion which governs the spatiai dependence of the acous- 
tic field for harmonic oscillations 

+ -0 , (1) 
where tp is the acoustic potential and k is the wavenum- 
berj the standard integral representation of the exterior 
potential is found to be^'® 

LI 

where the term 9/on, represents an outward normal 
derivative with respect to the body S as shown in Fig. 1, 
that is 


8<P(9) 

On, 


-V, (piQ)-n^ 


(3) 


Also, G{Pf (?) is a fundamental three-dimensional solu- 


From Eq, (2), If the acoustic potential and the normal 
acoustic velocity 9 ^(G)/9n, are known at each point on 
the surface of the body, then the acoustic potential may 
be calculated anywhere in the exterior domain . 

To solve Tor the surface potential, the point P Is 
moved to the surface of the body, and Eq* (2) then be- 
comes 

/./(*'<« - «'■'«> *”«■>'’> • 

(5) 

For the inhomogeneous Robin boundary Condition em- 
ployed In this study, a relatlbn between 9(fi{Q)/0n, and 
y?(Q) exists and is given by 

^-T(<?)?>(<?)=^A(0) . (6) 

so that Eq. (5) can be written in terms Of the potential 
only, that Is, 

^2n<pm+J^fA{!Q)G(J>,Q)d (7) 

If the acoustic Velocity A(Q) and the admittance F(Q) are 
specified at each point on the surface of the body, then 
the acoustic polential may be calculated at each point 
using Eq. (7). 

As mentioned earlier this equation does not yield 
unique solutions when the wavenumber k is an Internal 
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■2"(«>('’)*“^). 

where Bcp/Bti and <f> are related by Eq, (6). Equation (9) 
will yield unique solutions if the complex coupling con* 
slant is properly chosen. It is shown that a must meet 
the following restrictions to guarantee that Eq. (9) yields 
unique solutions^! 


m 


lm(of) #0, real or imaginary , 
Im(a) * 0 , k complex , 


(10) 


A problem arises in the numerical solution of Eq . (9) 
as the third term on the left hand side is strongly slngu* 
lar in its present form as the point Q approaches the 
point P on the Burface of the body. The authors of this 
paper have shown that this difficulty can be overcome 
by a proper interpretation of this singular term,' Em* 
ploying a vector transformation* and taking the Cauchy 
Principle Value, Eq. (9) is shown to be equivalent to 


clgonvaluo associated with the problem under considera- 
tion , Since these eigenvalues arc not known n p riori 
for general bodies, the formuiatlon cannot l>e relied upon 
to give consistently good results . 'n»ere are a number 
of papers in the ilteratnre*"Vdenling with this problem, 
and the relative mcrlle and shortcomings of the methods 
employed are discussed in detail in Itef. 1. 

An attractive approach from an analytical point of 
view is provided by Burton and Miller , * wlio liavc sug- 
gested the use of the following Identity to derive an al- 
ternative integral equation for the acoustic potential at 
the surface,* 

a, aa .//(yw) mm . . 

9M> JfJ V ®w« / 

( 8 ) 

This equation can now be solved for <p{P) by using Eq, 

(6) to relate the normal acoustic velocity and the poten- 
tial .at the surface, However this integral equation has 
its own set of associated eigenvalues at which unique 
solutions cannot be obtained. To circumvent the prob- 
lem associated with the solution of the integral equations 
derived from EqS . (5) and (8), Burton and Miller sug- 
gested the solution of the following linear combination 
of these equations : 

. , ; : . . . - : ; ,1 

fj{m -oiP.Q) (,W) -WP) 1 2 ^ rfs, 

-C"P(P)jJ fa, ■ »,)(«■)* C(P, QMS, -a JJ , 

All of the terms in Eq, (11) are now well defined; however, all Integrands are osctllatory and singular so tiiat 
care must be taken in their numerical approximation, 

B. Axisymmstric formulation 

When dealing with a body of revolution as Shown in Fig, 2 an axisyrnmetrlc formulation of the problem is advan- 
tageous « ‘° This being the case an element of area dS, becomes p ds dO, where s is the distance along the perimeter 
of the surface in the p -r plane , Assuming an .acoustic velocity distribution of the form 

~-su(s)co8(mS) , (12) 

and describing the $ dependence of the potential function by 

*[s)nip/co8i,m9) , ( 13 ) 

and letting S, = 0 (so that cos9,»l) Eq. (11) becomes: 

[♦(*«) cos(w9,) - +(s,)] r(s,)G(P, $)cos(wd,)dS, 

-o ff t/(s,) cos(m9,)dS,=2r{*(s,) + atr(s,)) . (14) 


( 11 ) 


Now, three sets of functloniS are defined: 

Influence functions: 

rc(p,Q) C08(we,)rfe, , 

•'0 


(15) 


Kernel functions i 
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Ki{r„)*i f co8(wg^)rfg, , 

* , ; * (16) 
^»(>».) * 2“ ■ cos(w9,) fl, 9, } 

Forcing functions: 


r,(f„).2o f G(P,(i)mH», • 




a?) 


where r,^ is tlio distance between points P and Q and 
and ttf are the outward normals to the surface at points 
i* and <?, respectively, In evaluating ff; and F^, the 
point at which is excluded from the Integration, 
Substituting Eqs, (15)-(17) Into Eq. (U) gives 

/ * - ♦(s,) f tF|(>V,) 

a •'O > 

vg 'A "■ 


where I is the length of the generating line of the surface 
of revolution, The $-(? coordiiuite directions have now 
heen essentially uncoupled so that the problem lias been 
reduced to the evaluation of the line integrals in the co- 
ordinate directions on the surface of the body. Tills 
formulation docs not restrict the form or type of 
boundary conditions on the body; it merely assumes that 
the boundary conditions can be represented by a sum 
(expanded in a set) of tangential modes. 


II, results 

The acoustic fields for a spheroi cylinder , and Inlet 
configuration have been computed by numerical solution 
of Eq, (18) using the tecfmlques described In Ref, 11, 
B«isicaUy, this method Consists of first specifying the 
p— z coordlnalos and the normal vector at each point on 
the surface. From these quantities the distances and 
the normal deriyatives can be obtained. The integral in 
Eq , (18) Is then separated into n integrals taken over 
subintervals of length l/n. The acoustic potential is as- 
sumed constant over each subinterval and the integra- 
tions are performed numerically using Gauss-Legendre 
quadrature In the p - z plane, Over the subinterval con- 
taining the pointPy the Integrand In Eq. (18) becomes 
infinite since approaches zeror Thus, only an even 
number of points is used in the quadrature algorlthlm, 



FIG* 2. Cylindrical surface geometry, 
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FiO* 3. Effect of the coupling Constant on the computed sur- 
face potential for a sphere of unit radius with 20 sublnlervalst 


since an odd number would necessitate inclusion of the 
point where A Gauss-Legendre quadrature for- 

mula is used in the circumferential direction to evaluate 
Eqs, (15)-(17). All calculations were performed on the 
Georgia Tech CDC Cyber 70/74 With 16 significant fig- 
ures. 


In all geometries investlgated| exact solutions were 
obtained for by assuming a monopplc source located 
at point (pjz) M (0^0) inside the surface. The normal ve« 
locitics and/or admittance values are then computed at 
each point on the surface using Eq, (6) and taken as the 
boundary conditions in Eq. (18), The surface potential 
4>(s^) is then computed from Eq, (18) and the farftold 
potential is obtained by numerically solving Eq, (21 with 
Eq, (6). The computed surface and farfield potentials 
are then compared with the known potential distribution 
of the monopole source 


<p(P)^-c^*«A , (19) 

where a is the distance from the source to the observa- 
tion point . For »l«l a dipole source was used to gen- 
erate exact solutions , and for m « 2 a quadrapoie source 
was used, 


To Investigate the effect of the coupling constant ot in 
Eqs, (15)-(17), the surface potential dlstrlbtiUons for a 
sphere of unit radius with a uniformly vibrating surface 
e , , m a 0) were computed for ot ^ 0, i, and i/k. Twenty 
subintervals were taken in the p -z plane, a four-point 
Gauss -Legendre quadrature formula was used over each 
subinterval and a 20-point Gauss-Legendre formula was 
used in the 6 direction , The magnitude of the potential 
should be unity at all points on the surface . The results 
presented in Fig . 3 show the computed magnitudeB of the 
surface acoustic potential to be in error by 12% for 
asi O at nondlmenslonal wavenumbers ka close to r, 2n, 
and 3r. These results are those that Would be obtained 


i 


from Eq , (7) . The relatively large errors arc expected 
from the analysis of Burton and Miller* and from previ- 
ous investigations using Eq. (7). Burton proves that 
setting the imaginary part of ot nonzero guarantees 
unique solutions to Eq, (18), For as? t the maximum 
error is reduced to less than 4% except when Af is close 
to 8, 0 > However, when and for sufficiently high 
values of ka, Eq. (9) is dominated by terms arising from 
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TABI^E It Effect of llio coupling pariimctcr a on the computed 
vuluoH of the nurfaco potential for a sphere. On the surface 
Y{0"O, V’M«("“i***»con8tnnt, w»0. Alt 
values of An correspond to internat clgcnfre<|ucndcs. 'Fwenty 
subintorvuls were taken In the p-c plane. 




0 

i/k 

i 

Exact 

n 

1>r 

2,6 

irboo 

0.0D8 

" 'l ” 


<i>i 

-0.3 

0.001 

-0,012 

0 

4,493409 


0. 190 

0.217 

0,208 

0,217 


n 

0.979 

0,970. 

0,956 

0,970 



-2,0 

-1.000 

-0.006 

-1 


Vt 

12,0 

0.000 

0.031 

0 

7.726252 

<Pr 

-0,081 

-0. 128 

-0.400 

-0.128 


n 

- 0, 994 

-0.992 

-0,872 

-0,992 

3jt 


2,0 

1.000 

0.005 

1 


Vl 

-19.0 

O.OOO 

-0,050 

0 


Eq, (8), As a result, the solution equations become ill« 
conditioned when ka Is sufficiently high and close to one 
of the elgcnfrequcnclos associated with the integral equa- 
tion based on Eq. (8), in Table I computed values close 
to these olgenfrcquencles and the cigcnfrequencies of Eq, 
(7) are compared with exact results for orsQ, I, and i/k, 
In all case, the value of i/k gives the most accurate re- 
sults, In Table II, the effect of introducing an admit- 
tance condition is presented for a » i/k, The admittance 
Y{Q) and forcing function A (<?) in Eq, (6) are chosen so 
that the relations 

| . ( 20 ) 

are eaiiafied on th^ surface and the exact solutions can 
be readily computed, Ttie loss in accuracy when an ad^ 
mlttancc condition Is used is minimal and restricted to 
the third significant figure. However, for discontinuous 
boundary conditions, where the forcing function is spe** 



FIG # 4i DependCncQ of the computed surface potenti^ for a 
finite cylinder with a zero admittance and nonzero normal Vclo-- 
city ever^'where on the surface for 20 subintervals. 
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FIG. 5< Effect of dtBCpntlnUOUB boundary conditions on the 
accuracy of the computed surface potential for a cylinder. 


cified over one part of the surface (l.e. , the admittance 
is zero there) and the admittance is specified over the 
remaining surfacoi errors of over ten percent in the 
real and Imaginary parts of the computed surface poten-^ 
tial result. For comparison , the case of a constant 
forcing function and admittance over the sphere for u 
is also presented and in all cases yields results of less 
accuracy than those obtained with a 

In this study consistently good results wore obtained 
with In Fig. 3 the computed and exact values 

for as; i/k agree to three significant figures over the 
range of nondimensiotu&i wavenumbers from one to ten, 
In fact, for this value of pt, the accuracy is significantly 
better at all wavenumbers investigated. While Burton 
and Miller* provide no reoommendattons for choosing 
one value of a over any other value with an Imaginary 
component, the choice used in the present study 

can be explained as follows . The terms In Eqs , (15)-* 



FIG . 6. Computed surface potential for a cylinder at the fl rat 
and second tangential modes for ka « 2. 
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FIG. 7> CumputccI farfiold potontlal distribution for n cylin- 
der ut m*0, and 20 radii from the center. 


(17) which involve a are of order k*, whereas the re- 
maining terms arc of order A*. Therefore, at higher 
wavenumbers the terms of order A* dominate. By choos- 
ing a to vary inversely With the wavenumber, all terms 
In Eqs. (15)-(17) remain of the same order with respect 
to wavenumber.’* 

A problem of more practical importance Is the finite 
axisymmotric duct since this surface approximates an 
engine configuration. The surface imtcntlal distribu- 
tions arc presented In Fig, 4 at different nondiniensional 
wavenumbers for wi bO, The normal acoustic velocity 
distribution A{Q) is chosen so that the solution for the 
acoustic potential satisfies Eq, (19), The parameter a 
is taken to be i/k. Twenty subintervals are taken in the 
p ^ z piunc and a 20-point Gauss-Legendre quadrature 
is used in the 0 direction. In Fig. 4 the variations of 
the magnitude and phase with distance along the perime- 
ter s are presented. The largest errors In the computed 


magnitude of the potential of about lty% occur on the ends 
of the cylinder and at the corners, The results at the 
ends can be improved without Increasing the number of 
points by area weighting rather than taking equidistant 
points along the perimeter, The errors in the phase are 
less than 4% in all cases. The errors in the magnitude 
of the computed surface potential Increase with increas- 
ing nondlmcnsional wavenumber, but, even when A'(t> 10, 
the numerical results are within 10% of the exact solu- 
tions, For a «o or f the errors are significantly larger 
above A(i«2. 

» 

In most inlet problems the boundary conditions are 
discontinuous with the acoustic, velocity or potential 
(which Is directly pro|x>rtional to the acoustic pressure) 
specified over part of the surface and the admittance 
frepresenting liners) over the rest. To determine the 
o;fcct of the discontinuities and the use of an admittance 
function on the numerical results for a cylinder 
was investigated, The velocity was specified on the 
ends and the admittance was specified in the center so 
that the soiction for 4> was given by Eq, (10) and Eq, 

(6) is satisfied. Again, 20 points are used In the p-s 
and B directions. The results are shown in Fig. 5. Al- 
though the errors in the numerical results for this case 
are higher than those observed in Fig. 4, the errors 
still remain within 10% for values of k<t less than 5. 
However, when ka * 10 errors of up to 40% in the magni- 
tude of the potential are encountered close to the dis- 
continuity in the boundary condition. This error can be 
reduced by Increasing the number of subintervals In the 
p - 2 plane. Doubling the number of subintervals halves 
the error, When both the normal acoustic velocity and 
the admittance arc continuous on the surface, the errors 
are of the same order of magnitude as those of Fig, 4 . 
For tangential modes, the variation In the circumferen- 
tial direction behaves as cosmS where 1,2, , , , . To 
check the numerical integration scheme in the circum- 
ferential direction, the surface acoustic potential was 
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computed for mi rI and hi *2 for the cylinder shown In 
Fig. 4 4 Tfic results are presented in Pig. 6 for ka»Z 
with the normal acoustic velocity spocincd everywhere 
on the surface, The computed and exact results (I ,e, , 
from a dipole and quadraiiole) arc In jrsrccment to within 
2% for both HI and HI ° 2, 

It has been shown* that once the surface potential has 
been accurately computed, the farfleld can bo deter- 
mined to at least the same accuracy as the surface po- 
tential, This result Is confirmed by the data presented 
In Fig, 7 for the cylinder of Fig, 4 with the velocity 
specified everywhere on the surface With A*a « 2 and ni * Q. 
The results at 20 radii from the surface arc in agree- 
ment with exact results obtained from Eq, (10) to within 
1% even though the surface errors at some points arc 
above 2%, Data In Fig, 8 show that accurate results are 
obtained at distances greater than one integration step- 
size from the surface. At closer distances errors from 
the numerical evaluation of the singularity In the Green’s 
function defined by Eq, (4) leads to large errors. 

The studies of the acoustic fields of the sphere and 
cylinder Served to check out and refine the numerical 
procedures and programming techniques . The next con- 
figuration Investigated was an inlet used in a study l>y 
NASA*” 'HUs Inlet Is shown In Fig. 8 and was chosen 
because; 

(1) unlike most inlets used in research studies, it 
docs not have a bctl-mouth shape but is shaped like a 
typical Inlet used In existing aircraft, and 



FIG . 10. Eltept of Inoreaslng the number of aubintervals in 
computtng the surface potential lor the Inlet conffguratlon at 
Jfo R 1, m ■= 0, 
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(2) complete details on generating the Inlet boundary 
are given In Ref. 13. For this Inlet, all cases were in- 
vcstlgated with a»t/k, 

As Btcn in Pig. 10, the normal velocity distribution, 
which reprosents a forcing function, Is highly dlscon- 
tinuoUB and provides a severe test of the numerical tech- 
niques employed. The numerical and exact solutions 
for the surface acoustic potential are compared in Fig. 

10 for 32 and S4 sublnlervals taken along the perimeter 
of the inlet In the p z plane , Because of the errors in 
approximating the lengths of each subintcrval, the exact 
solutions differ slightly as the distance along the perime- 
ter s increases. The contorbody in Fig. 0 extends from 
0^ $£0*6, the fan inlet covers 0.64 1.4, the interior 
contour extends from 1. 4 4 3. 6, the exterior from 

3,54 s £5, 5, and the circular arc lies within the inter- 
val 5.54s47,45, Increasing the number of points de- 
creases the error proportionately as indicated by the 
data in Fig. 10 at a nondimcnslonal wavenumber ka of 
unity, where a Is the radius of the inlet at the fan entrance 
section. The absolute average error in the results de- 
creases from 10,2% for 32 subintervals to 4> 16% for 53 
sublnlervals. The computation time increased from 53 
to 143 s, respectively. 


TABLE II. Effect of specifying nn admittiuicp on the computed 
surface potential for a sphere. In all cases twenty sub- 
lnlervals arc taken in the p -x plane, and every-, jjl 

where on the surface. For Case I, *4 (<?)•><’'* and I'W) •’t^wery- 
whore on the surface. For case IQ, and FW) 

»0 over 1/8 of the gurface andAW)* 0, Y^(l-ik) over the 
remainder, Caac Q la considered in Table 1, 


— — 

cat 

lel 

Cuflel 

Caiie 11 

Case m 

Exact 

ha 

a* 

'A 

o#0 

«*f/* 


Values 

1 


-0.539 

-0.637 

-0.638 

-0.52 

-0.540 


n 

-0.845 

-0.849 

-0.843 

-0.87 

-0,842 

2 


0.416 

0.422 

0.417 

' 0.43 

0.416 



-0*911 

-0.937 

-0*909 

-0.92 

-0.909 

2 


0. 093 

0,916 

0.990 

1.00 

0.990 


n 

-0*142 

-0.496 

-0.140 

-0.16 

-0,141 

5 


- 0. 285 

-0.288 

-0.284 

-0.25 

-0.284 


n 

0. 961 

1. 145 

0.959 

1.00 

0,959 

10 


0. 641 

-0.3 

0.639 

0.90 

0,839 



0.540 

0,9 

0.544 

0.49 

0.544 
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FKJ, 1 1 . Effect of tnorenstiK frequoney for tho inlet Configura- 
tion at Mt « 0 , n « 53 on the comiMtod aurfneo potential. 


As siiown In Fig. 11, the crrorB incrciiie With In- 
creasing frequency. Like the cylinder, the maximum 
error in the acoustic potential for the inlet configuration 
occurs at the points of discontinuity, The average error 
inercasoB from 4.16% at ka^l to 15% at A’fl*10, 

For the data in Figs, to and It, the acoustic potential 
is assumed constant in the tangential plane. The results 
for a cos(m6) distribution are presented in fig, 12 at 
ka^2, These results show the insensitivity of the accu- 
racy of the computed results to tho tangential distribu- 
tion for m * 1,2. The exact solutions were again gener- 
ated by assuming dipole and guadrapolc sources located 
at (p,i)»(0,0). 

Based on the above reauUaour numerical and program*^ 
mlng techniques (tre capable of yielding rell^le results for 
arbitrary gcpmctrles and boundary conditions* At higher 
frequenclesi It appears that more points must be 

ttiken to Increase the accuracy of the computed results. 

tii. SUMMARY AND CONCLUSIONS 

An Integral solution of the HGlmhoUz equation is 
veloped for use In acoustic radiation problems . Unlike 
the classical formulation which can lead to integral equa- 
tions that do not have unique solutions at frequencies 
corresponding to certain Internal eigenfrequencles of the 
region enclosed by the surface under consideration^ the 
formulation used in this study is valid at all irequencles • 
Also I unlike most cu r rent methods and for mulations it 
is stralghtfoinvard to implement regardless of how com- 
plicated the surface or the boundary conditions may be. 
The surface potentials computed numerically for a sphere 
and cylinder using 20 sublntervals along the perimeter 
and for an inlet configuration with 53 subintervals are 
accurate to within 10% for nondimenslonal wavenumbers 
ka of from one to ten, where k Is the wavenumber and a 
is the characteristic length. For discontinuous boundary 
conditions, the numerical and o?cact values are in agree- 
ment to within 10% for < At higher frequencies the 
results are as much as 40% in error at points of discon- 
tinuity, which suggests taking more points in evaluating 
the integral equation to increase the accuracy when dis- 
continuous boundary conditions are specified. Increas- 
ing the number of subintervals decreases the error pro- 
portionately. At distances greater than the numerical 
integration stepslze, the farfiold results are at toast as 
accurate as the corresponding surface potential solutions, 
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FIG , 12, Effect of mode number m on the computed Murtace 
potenUai of the Inlet configuration for Ae »2 and 53. 
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Fillwl for drawing our attention to the following paper* 
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